Effect of magnetic field on impurity bound states in high-T c superconductors 
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We consider the influence of a magnetic field H on the quasiparticle bound states near scalar 
impurities in d-wave superconductors. A "Doppler shift" in the excitation energies induced by the 
supercurrent leads to several important effects. At large but finite impurity strength, there are 
corrections to the energy and width of the impurity-induced resonance, proportional to H 2 . On the 
other hand, in the limit of very strong impurity potential (unitary limit), the bound state is destroyed 
and acquires a finite width proportional to H(\aH)~ . There are also considerable changes in the 
asymptotic behaviour of the bound state wave functions. 

PACS numbers: 74.72.-h, 74.20.Rp 



High-temperature superconductors (HTSC) belong to 
a class of unconventional superconductors with d-wave 
symmetry [n| . A non-trivial orbital structure of the order 
parameter, in particular the presence of the gap nodes, 
leads to the effects of disorder in HTSC being much 
richer than in conventional materials. For instance, in 
contrast to the s-wave case, the Anderson theorem does 
not work and non-magnetic impurities exhibit a strong 
pair-breaking effect. Also, a finite concentration of disor- 
der produces a non-zero density of quasiparticle states at 
zero energy, which results in a considerable modification 
of the thermodynamic and transport properties at low 
temperatures J2J. 

One of the striking features of d-wave superconductors 
is that even a single scalar impurity has a notable effect 
on the superconducting state, creating a quasiparticle im- 
purity bound state (IBS) in its vicinity whose energy and 
width tend to zero in the limit of strong impurity poten- 
tial H . Various aspects of the theory of IBS have been 
elaborated by many authors, both analytically P-p| and 
numerically ||. From the practical point of view, the 
high sensitivity of IBS to the symmetry of the order pa- 
rameter and the impurity strength makes these states a 
powerful tool for probing the properties of HTSC. This 
has been done recently in beautiful scanning tunneling 
microscopy (STM) experiments on BSCCO compounds 
|0-^|, in which many of the theoretically predicted fea- 
tures of IBS were observed, such as (i) sharp peaks in the 
energy dependence of the single-particle density of states 
(DoS), (ii) a cross-shape anisotropy of the quasiparticle 
wave functions in a-b plane, and (iii) a slow power-law de- 
cay of the wave functions away from the impurity. (Some 
apparent discrepancies between the theory and the exper- 
iment, such as the position of the maximum local DoS 
directly above the impurity site as seen in STM pictures, 
can be attributed to the blocking effect of Bi-0 layers 
| fiO| and thus do not invalidate the standard theoretical 
picture of IBS.) 

All previous studies have neglected the influence of an 
external magnetic field on IBS. It is known, however, 
that in the presence of the gap nodes, the supercur- 
rent induced by a magnetic field acts as an effective pair 



breaker and gives rise to a finite density of bulk quasipar- 
ticles at low energies. This should affect the properties 
of IBS because of the hybridization of localized and de- 
localized states. The purpose of the present article is to 
study the combined effect of magnetic field and strong 
impurity scattering in d-wave superconductors. We ne- 
glect the Zeeman splitting and concentrate on the phys- 
ical properties relevant to the c-axis STM experiments 
which probe the local DoS at various distances from an 
impurity site. The influence of the magnetic field on 
the disorder-averaged density of states in d-wave super- 
conductors with randomly distributed unitary impurities 
was studied in Rcf. ]ll| ]. 

Let us consider a repulsive point-like impurity with po- 
tential U(r) — ud(r) (u > 0) in a two-dimensional d-wave 
superconductor. The external magnetic field can be di- 
rected either in the a-b plane or along the c axis. In both 
cases, we assume that the supercurrent p s is locally uni- 
form ; and use the gauge in which the order parameter is 
real |l2| . The model Hamiltonian can be written in the 
following form: 



H = J2 c ^ 
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where Ck = (ck|,c_ k ,) T are Nambu operators, and 
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is the Bogoliubov-de Gennes (BdG) operator. The nor- 
mal part of -ffkk' depends on the supercurrent via the 
gauge transformation: H± = £k±p s + E/k-k' ~ £k i 
vj?p s + C/k-k', where £k = k 2 /2m — /j, is the normal 
state spectrum (we assume a spherical Fermi surface, 
which does not restrict the generality of our results, but 
considerably simplifies the calculations), Vf = Vk£k is 
the Fermi velocity, and C/k is the Fourier transform of 
the impurity potential. The mean-field order param- 
eter, corresponding to d x 2_ y 2 symmetry, has the form 
Ak = Ao(/c^ — hy) — Aq cos 2ip, where ip is the azimuthal 



angle in a-b plane. Note that we do not calculate the or- 
der parameter self-consistently and assume Ao to be con- 
stant. The numerical investigation of the self-consistent 
BdG equations shows that there is some suppression of 
the order parameter magnitude near an impurity site [pL3|| , 
which can affect the low-energy behaviour of DoS in a 
system with a finite number of impurities H] . Here we 
study the single- impurity limit, which corresponds to a 
dilute disorder concentration. In this case, it has been 
demonstrated in Ref. p5[ that the order parameter vari- 
ation leads only to renormalization of the effective impu- 
rity strength towards the unitary limit. 

The quantity measured in STM experiments is the lo- 
cal differential conductance dl/dV, which is proportional 
to the local density of states: 



N(r,u>) = --ImGf^rjw) 

7T 



(3) 



where G R is the retarded Gor'kov-Nambu matrix Green's 
function. In the presence of a single scalar impurity, one 
can express G R in terms of the Green's function of a clean 
superconductor. Thus, 

G R (r U T 2 ;u>) = G R ( ri -r 2 ,uj) 

+G R (r l ,u>)T(uj)G R (-r 2 ,u>), (4) 

where the T-matrix is given by T(w) = 
UT 3 [l-ug(u>)T3]" , and g(w) = G r (0,uj) = 
Jd 2 k/(2Tr) 2 G R (k,oj). The Green's function G R de- 
scribes a homogeneous system without an impurity, but 
in the presence of a uniform supercurrent. The effect 
of the supercurrent in a translationally invariant system 
amounts to a "Doppler shift" in the quasiparticle energy. 
In the momentum representation, we have 

(ui+ - v F p s )r + £kT 3 + A k ri 



G£(k,u) 



(5) 



( W+ -Vjrp.) 2 -£-A2 ' 

where u+ = lu + iO, and tj are the Pauli matrices. 

The energies of the impurity-induced IBS correspond 
to the poles of the T-matrix and satisfy the equation 
detT(w) = 0, whose solutions can be complex. We are 
interested in the limit of small supercurrent and strong 
impurity scattering, so that the relevant energies are ex- 
pected to be small compared to the gap magnitude. This 
case is of particular interest because the most profound 
effects related to IBS have been observed in the vicin- 
ity of Zn impurities in BSCCO, which have the s-wave 
phase shift Sq ~ 0.487T and thus are very close to the 
unitary limit M . The momentum integrals can be easily 
calculated, giving the following result at Imoj = 0: 



Img(oj,p a ) =g (uj,p s )T Q + gi(w,p a )n, 
where, in leading order in (vfPs/Aq) 2 , 
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In these expressions, Np is the DoS in the normal state 
at the Fermi level |[(|, i labels the gap nodes, and Vj are 
the Fermi velocities at the nodes. We see that, at small 
supercurrents, it is possible to neglect the off-diagonal 
terms in (|6|) (retaining these terms would give the cor- 
rections of the order of (vfPs/^o) 4 to the results below). 
From the Kramers-Kronig relations, we obtain 



Re#H = -— f- V (w-v i p a )ln 1 ^— -r tq. (8) 

2A f-f |u/-v,p s | 



The next step in the derivation of the T-matrix is to 
continue g(u>), whose real and imaginary parts at the real 
axis are given by Eqs. (||) and (||), to the whole complex 
plane of u. Introducing the notation z = u)/Aq, we have 

g(u>) = ttNfF(z)t , where 



F{z) = ^-Y,{z- Zl )Hz- Zl )-iz (9) 

1=1 

with Zi — Vip s /A (J2i z i — 0)- Four logarithmic branch 
cuts go down from z = z% parallel to the negative imagi- 
a.xis_ Fina.llv 



nary axis. Finally, 
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where c = (iruNp ) 1 = cot 5o > for a repulsive im- 
purity. This expression for the T-matrix is valid at 
\z\,\zi\ <1. 

The IBS spectrum is determined by the equation 



F{z) = ±c. 



(11) 



In the absence of a supercurrent, F(z) — ► Fq(z) — 
(2/tt)z In z — iz, and Eq. (O) can be easily solved at c <C 
1, giving the IBS energy of the form luq = zoAo, where 
Rezo = = F 7r c/2|lnc|, Imzo = — 7r 2 c/41n c with logarith- 
mic accuracy |3| . The presence of a non-zero imaginary 
part indicates that the impurity-induced bound state is 
in fact a narrow resonance. 

At p s ^ 0, the effect of supercurrent on IBS strongly 
depends on the relation between the Doppler shift VfPs 
and the "bare" energy ojq. If p s is parallel to one of the 
crystallographic axes, e.g. (100), then zx,4 = — z i,z — 
z s = vfPs/V^Aq, and Eq. (||) takes the form 
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n 



s=l 



F(z) = —iz -\ (z + z s ) ln(z + z s ) 

IT 

H — (z- z s )ln(z- z s ). 
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(12) 



For z s <C |^o 1 7 the solution of Eq. (|ll| ) can be found per- 
turbatively in z s . Using the expansion F(z) — Fq(z) + 
z 2 /irz, we find 



Re z = Re z =p — z^ 



Im 2 = Im Zq 



1 



2c In 



-zl. 



(13) 



These expressions are valid as long as \5z(z s )/zo\ <C 1. 
Thus, in the presence of a small magnetic field, the cor- 
rections to the IBS energy and width are proportional to 
H 2 (z 2 s = 0.5[H/H C (0)} 2 for H _L z). 

More interesting is the opposite limit of "large" super- 
currents, which is relevant when the impurity scattering 
is close to the unitary limit. In zero field, the bare IBS 
energy vanishes at c = 0, so the impurity-induced reso- 
nance gets decoupled from the continuum of propagating 
excitations and becomes a true bound state with zero 
width. The dominant energy scale in this case is pro- 
vided by the Doppler shift vfPs , which makes it possible 
to treat the right-hand side of Eq. ( |ll| ) as a small per- 
turbation. It can be checked that the equation F(z) — 
has only one solution in the complex plane: 

(14) 



2|lnz fl | 

At c 7^ 0, we look for a solution in the form z = z* +5z(c) 
and obtain 



Rez = =F 



Imz = 



2|lnz s 
2|lnz, 



+ 0(c 2 ). 



(15) 



These expressions are valid as long as | Rez/z»| <C 1, i.e. 
at z s 3> c. From Eqs. (15) we see that the zero-energy 



IBS in the unitary limit is destroyed by magnetic field, 
getting replaced by a resonance whose width is propor- 
tional to if (In H)^ 1 . The physical reason for this is clear 
from Eqs. (|3j) and (J6J): at non-zero supercurrent, the 
DoS of the bulk excitations does not vanish at u> = 0, 
which leads to a stronger hybridization of the IBS and 
the continuum of pro pag ating states. We would like to 
stress that the result ( |l5[ ) is manifestly non-perturbative 
in magnetic field. 

The calculations for other directions of supercurrent 
can be done in a similar fashion. Solution of Eq. (fnl) 
leads to the expressions analogous to ([13]) and ([l|), al- 
beit with different numerical coefficients. Therefore, the 
qualitative effect of magnetic field on the IBS energy does 
not depend on the direction of p s . 

In order to visualize our results and facilitate the com- 
parison with STM experiments, we have computed the 
local DoS at the impurity site: 



N(0,u>) = -JVpIm 



cF(z) 
■-F{z) 



(16) 



where F is given by Eq. (||). The dependence of N(0,u>) 
on energy and supercurrent is plotted in Fig. [y. We see 



that the asymmetric peak in the local DoS, which corre- 
sponds to a hole-like resonance at c > 0, gets shifted and 
broadened in the presence of supercurrent. 

It follows from (|l|) that N(0,u>) = at c = 0. For 
this reason, the STM measurements directly at the im- 
purity site are not useful in the unitary limit. To study 
the effect of supercurrent on IBS in this case, one should 
calculate the local DoS away from the point r = 0, e.g. at 
one of the nearest neighbors of the impurity site, where 
the local DoS is known to reach its maximum value at 
p s = @. According to Eqs. (§), (|) and ©, the local 
DoS at site r in the unitary limit can be represented in 
the form 



N(t,w) 



1 



hn-—[G^(r,cj)G^(-r,u;)} ll , (17) 



TT 2 N F F{z 



where Nq(u>) = —(l/n)go(oj) is the DoS for a clean d- 
wave superconductor in the presence of supercurrent, 
and Gq(t,uj) is the Fourier transform of Eq. (g). At 
r = a and (uj,vfPs)/^o —> 0, one can neglect the first 
term on the right hand side of Eq. (|17| ) compared to 
the second one, which is singular in this limit. The sin- 
gularity comes from F^ 1 (z), while the product of two 
Green's functions is not singular and can be replaced 
by its value at u = p s = 0, which is real. There- 
fore, N(sl,uj)/N f ~ ImF _1 (z). We have calculated 
the energy dependence of 8N(a.,u) — N(a,uA -~ No(oj) 
for p s || a and plotted the results in Fig. g. As the 
supercurrent increases, so does the width of the zero- 
energy peak, whereas its magnitude decreases. Exper- 
imentally, a notable suppression of the c-axis zero-bias 
conductance peaks (ZBCP) has been reported in Ref. 
for YBCO/Ag and TBCCO/Au planar junctions. It 
is tempting to attribute this observation to the effect of 
magnetic field on the IBS induced by strong defects at 
the surface of those materials. 

Another peculiar property of IBS, which was predicted 
in Ref. || and observed in recent STM experiments J7[g[ , 
is a sharp four-fold anisotropy of the IBS wave functions 
at zero bias, with a characteristic l/r 2 -tails at large dis- 
tances from the impurity. This not only reflects the mi- 
croscopic symmetry of superconducting state, but also 
has important consequences for the quasiparticle trans- 
port in HTSC. It was argued in Ref. |18| that, in the 
presence of many impurities, the overlap of the extended 
IBS wave functions along the gap node directions can re- 
sult in the formation of an impurity band, where all states 
are delocalized and can thus participate in quasiparticle 
transport. This should be contrasted to the propagating 
excitations which are localized in the presence of short- 
range disorder potential |l!j . 

Here we address the following question: What is the 
effect of magnetic field on the asymptotics of the IBS 



wave functions in the unitary limit? We study the de- 
cay of SN(t,cj) — N(r,oj) — N (uj) along the nodal di- 
rection (110) at w = 0. The coordinate dependence of 
Gjf(r, u>) in Eq. ( |l7| ) is described by rather cumbersome 
expressions (the details of calculations will be given in a 
separate publication), which take on a particularly sim- 
ple form in the limit of large distances from the impurity 
r 3> pj 1 S> £o — ^f/Ao- In contrast to our results for 
the IBS energies, the asymptotic behaviour of the wave 
functions strongly depends on the direction of the su- 
percurrent. If p s is perpendicular to the nodal direction 
(110), then, in leading order in 1/r, 



SN(r,0) 



j(Ps)N F £$ sin 2 k F r 



4tt 2 



for k F £,o > p s r > 1, and 



6N(r, 0) 



l{Ps)Ps 1 

" 8tt 2 A r' 



(18) 



(19) 



for p s r 3> k F t;o 3> 1 (we introduced the notation "f(p s ) 
ImF -1 ^ = 0)). If p s is parallel to (110), then 



6N(r,0) = - NfP '$ hn 



F-^O)- 



Jlih-pr ' 



(20) 



Eq. (jl8| ) is essentially the zero-field result, whereas Eqs. 
( |l9| , l20D show that magnetic field leads to a slower power- 
law decay of the envelope IBS wave functions along the 
nodal directions, or even the disappearance of the Fermi 
oscillations at very large distances from the impurity. 
These changes reflect the contributions from the nodal 
quasiparticles excited by supercurrent. In particular, the 
non-oscillatory tail in Eq. (|19| ) comes from the (110) and 
(110) nodal states for which ki?r = 0. 

In conclusion, we have studied the influence of mag- 
netic field on the impurity bound states in rf-wave su- 
perconductors. We have found several effects, which 
can be directly measured in STM experiments, such as 
a non-linear shift of the electron-hole asymmetric peak 
in the local DoS at c ^ 0, and also strong suppression 
and broadening of the zero-bias peak in the unitary limit. 
The changes in the wave function asymptotics can lead to 
a stronger long-range overlap between the bound states 
at different impurities, which might considerably affect 
the quasiparticle transport in HTSC. 
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FIG. 1. Evolution of the DoS at the impurity site as a func- 
tion of energy at increasing magnetic field (p s || a, c = 0.1). 
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FIG. 2. Suppression and broadening of the impu- 
rity-induced DoS peak at a nearest neighbor site in the unitary 
limit at increasing magnetic field for p 3 | a (a sharp peak at 
zero field is not shown). 



